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Adaptive perturbation is a new method for perturbatively computing the eigenvalues and eigen- 
states of quantum mechanical Hamiltonians that heretofore were not believed to be obtainable by 
such methods. The novel feature of adaptive perturbation theory is that it decomposes a given 
Hamiltonian, H, into an unperturbed part and a perturbation in a way which extracts the leading 
non-perturbative behavior of the problem exactly. This paper introduces the method in the context 
of the pure anharmonic oscillator and then goes on to apply it to the case of tunneling between both 
symmetric and asymmetric minima. It concludes with an introduction to the extension of these 
methods to the discussion of a quantum field theory. A more complete discussion of this issue will 
■ be given in the second paper in this series. This paper will show how to use the method of adaptive 

' perturbation theory to non-perturbatively extract the structure of mass, wavefunction and coupling 

C . constant renormalization. 

PACS numbers: 03.65-w,03.70+k,11.10 

CN ■ I. INTRODUCTION 

> 

ly-j ' Perturbation theory, one of the most useful tools in the physicist's toolbox, is the method most commonly used to 
analyze problems in quantum mechanics and field theory. It is so useful that it is applied even when the perturbation 
parameter isn't small enough for the expansion to converge. One then typically tries to extract meaningful results 
by applying some summation technique to the resulting series, often with mixed success. What is really needed is a 
1^ • method that extracts the leading non-perturbative behavior of the quantity to be computed, and then corrects this 
leading approximation in a systematic manner. This paper introduces a simple technique that works in just this way. 
I call this method adaptive perturbation theory 1\. 

I will limit the bulk of this paper to discussing how adaptive perturbation theory works on problems in quantum 
Oh, mechanics. In the second paper in this series I will show how to apply an appropriately generalized version of adaptive 
^ ' perturbation theory to a quantum field theory. In particular, I will then show how to extract, in a non-perturbative 
context, the general structure of mass, wavefunction and coupling constant renormalization. 
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II. EXAMPLE 1: THE PURE ANHARMONIC OSCILLATOR 

To introduce the ideas behind adaptive perturbation theory, I begin by studying the pure anharmonic oscillator 
since, naively, this problem seems to be impossible to treat perturbatively because it is devoid of a small expansion 
parameter. The pure anharmonic oscillator is defined by a Hamiltonian of the form 

H=^-p'+^-Xx\ (1) 

where the operators x and p satisfy the usual canonical commutation relation [x,p] = i. 

Since H contains only the terms p^ and A , there does not seem to be a way to perturbatively compute the 
eigenstates and eigenvalues of H based upon a division of H into an unperturbed Hamiltonian Hq, whose eigenstates 
are known, and a perturbation V. Nevertheless, I will now show that this can be done. The trick is to define a 
decomposition of the Hamiltonian in Eq. ^111^0 an unperturbed (i.e., solvable) part and a perturbation, in a way 
that doesn't depend upon A. The choice of decomposition is what is adapted to the problem at hand. Thus, the 
decomposition, and therefore the perturbation theory for each eigenstate, will be different. 
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To begin the process of decomposing H, I introduce a one-parameter family of annihilation and creation operators, 
Al^ and A^, as follows: 




Note that given the commutation relations of x and p, the commutation relation [A^,/!^] = 1 holds independent of 
7. In terms of these operators H can be rewritten as 




where I have defined the 7-dependent number operator to be 

= A\A^. (4) 
At this point I define a 7-dependent decomposition of H into an unperturbed part and a perturbation as follows: 

H = Ho{j) + V{j), 
Hoh) = (j-f^)(2iV, + l) + A^,(iV,-l), 

+ 4^ {liA" + A^) + + 4 A')) • (5) 

In other words, the unperturbed Hamiltonian is defined to be that part of H which is diagonal in the number operator 
Nj and the perturbation is defined to be everything else. 

If we think of the Hamiltonian as a large matrix in the basis of Fock-space states generated by applying Al^ to the 
state |0-y), then -ffo(7) is the diagonal part of the matrix and the perturbation is the matrix obtained by setting all 
terms on the diagonal to zero. 

Although it isn't conventional to include the term -^N^ {N^ — 1) in the definition of the unperturbed Hamiltonian, 
Ho (7), it makes sense to do so for two reasons. First, this part of the Hamiltonian is diagonal in and thus it is 
just as easy to deal with as the part proportional to the first power of Nj. Second, including it in the unperturbed 
Hamiltonian puts a term proportional to 7V^^ in the energy denominators which appear in the perturbation expansion 
based upon this decomposition. This suppresses the n! growth of the usual perturbation expansion. The only apparent 
drawback to putting a term proportional to N^'^ into Hq{^) is that it makes the energies of the unperturbed levels 
grow like n^, whereas the correct answer is n^^^ . While this is not a disaster, it seems to say that getting the correct 
answer from a straightforward perturbative approach will be difficult. This observation leads us to add the adaptive 
element to adaptive perturbation theory. 

The essential point is that I have not as yet committed to a particular value of the parameter 7 in Eqs. I5I5I The 
trick consists of adapting the choice of 7 to the computation at hand. Thus, to compute the energy of the N*-^ 
eigenstate I will choose a value, j{N), which varies as a function of N. The appropriate value of 7(A^) is fixed by a 
simple variational calculation. 

The variational calculation, which serves as the starting point for any adaptive perturbation theory computation, 
begins by defining a 7-dependent family of Fock-states, \Nj). First we define the 7-dependent vacuum state, |0^), by 
the condition 



A^\0^) = 0. 



(6) 
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Next, we define the n-particle state built on this vacuum state such that 

Nj\n^) = n|n-y), (7) 

by 

K) = 47<'|07>- (8) 

In terms of functions of x, the state |0^) is just a gaussian whose width is determined by 7. Similarly, the state 
is the appropriate n*'' order Hermite polynomial obtained by applying the 7-dependent differential operator 

to this gaussian. (It is worth noting that the specific functional form of the state depends upon 7, but the number 
of times this function takes the value zero does not.) 

The value of 7 used to define the adaptive perturbation theory for the n*'' level of the anharmonic oscillator is 
determined by requiring that it minimize the expectation value 

En{j) = {n^\H\n^). (10) 

Eqs. l5l5l show that this expectation value is equal to 

2 + ^) (2n+l) + A 
4 472/ ^ ' 472 



£;„(7) = (A^^|i?o(7)l«>= 7 + 7-2 (2^1+1) + T-2n(n-l). (11) 



Minimizing En^j) with respect to 7 leads to the equation 

At this point all we have to do is substitute this value into Eg. 1101 to see why we chose 7 to minimize i?„(7). If we 
do this, we obtain 

EnilUn = I (2n + 1)'/' (2n2 + 2(n + 1))'/' (13) 
o 

which, for large n, behaves as X^^^ n^^^ , which is the correct answer. 

The fact that all energies scale as A^/^ is an exact result and so the non-trivial part of variational computation is 
the derivation of the dependence of the energy on n. To see why all energies are proportional to A^/^, it suffices to 
make the following canonical transformation 



A^/V (14) 



Ai/6 

In terms of these operators, the Hamiltonian of the pure anharmonic oscillator becomes 

i/ = AV3Qp2^i,4^, (15) 

thus proving the claim. A comparison of the variational computation and the result of a second-order perturbation 
theory for A = 1 and widely differing values of N is given in Table As advertised, the adaptive perturbation theory 
for each level converges rapidly, independent of A and N. 



A. Why does the expansion converge? 



At this point I should address the question of why adaptive perturbation theory converges. 

First, observe that once I have done the variational calculation, the perturbation theory computation is done for 
states defined in terms of a single fixed value of 7. Thus, all states which appear in the perturbative expansion 
are orthogonal to one another and the perturbation series is well-defined. Second, for the values of 7 associated 



4 



with a given state \no), the unperturbed energies for states having another occupation number, n, all have a term 
« r? lr?J'^ . It is only for values of n that are near no that the energy denominators behave as n^^^ . For values 
of n >> no the factor 7 does not suppress the n^ behavior. This is what suppresses the n!-growth of the perturbation 
expansion, since now the factor of •\/(n + l)(n + 2)(n + 3)(n + 4), which comes from applying to a state |iV-,), is 
canceled by the 1/n^ appearing in the energy denominator, instead of the usual 1/n. I'll make no atternpt to prove 
this assertion, instead I refer the reader to a pair of very nice papers by I. G. Halliday and P. SuranyiU 01 where, 
by a very different technique which shares this same feature, they give a proof that the fact that the unperturbed 
energies are quadratic in iV leads to a convergent perturbation expansion. 

Finally, I should comment that the fact that variationally determined states |n^) are not orthogonal to one another 
isn't a problem. First, there is a simple formula that expresses the state |0-y) as a sum over all states |n-y') when 7' 
different from 7 (see Appendix A). Thus, if one wishes to do a non-perturbative diagonalization of the full Hamiltonian 
truncated to many different variationally determined states, one can either orthogonalize these states by hand, or use 
the formula for the overlap of these states and solve the so-called relative eigenvalue •problem 

H\i^)^tM\^). (16) 

Note, that if |n-y(„)) stands for one of the variationally constructed states, then the Hermitian matrix M is given by 

M„,™ = (A/^(„)|M^(„)). (17) 

The second point is that since the perturbation expansion for any state |n^(„)) converges, eventually the perturbatively 
corrected states |n^(„)) become orthogonal to one another. This follows from the fact that eigenstates of a Hermitian 
Hamiltonian having distinct eigenvalues are orthogonal to one another. 

B. What does this have to do with quasi- particles? 

An interesting corollary to the adaptive perturbation theory technique is that it provides an explicit realization of 
the quasi-particle picture underlying much of many-body theory. What we have shown is that no matter how large 
the underlying coupling, the physics of the states near a given no-particle state can be accurately described in terms 
of perturbatively coupled eigenstates of an appropriately chosen harmonic oscillator. Of course, as we have seen, 
the appropriately chosen harmonic oscillator picture changes as no changes. Furthermore, the perturbatively coupled 
states \n^{nQ))i for n ~ Uq-, correspond to states containing an infinite number of particles, if we choose as a basis 
those states that correspond to a significantly different value of no- 

III. EXAMPLE 2: ADDING A MASS TERM 

An obvious modification of the pure anharmonic oscillator is to add a positive mass term; i.e., let iJ be 

H =\v^ ^"^.^ (18) 

This problem is dealt with in exactly the same way. First, introduce 7-dependent annihilation and creation operators, 
to arrive at a slightly modified formula for the expectation value in state |iV-y); i.e., 

£;„(7) - (J + f + (2 n + 1) + - !)• (19) 



TABLE L A comparison of the zeroth order and second order perturbation results, for the energy of the A*'* level of the pure 
anharmonic oscillator, to the exact answer for A = 1 and widely varying values of A. 



A 


A 


Variational 


2"'* Order Perturbation 


Exact 


Variational % Err 


Perturbative % Error 


1.0 





0.375 


0.3712 


0.3676 


0.02 


0.0098 


1.0 


1 


1.334 


1.3195 


1.3173 


0.01 


0.0017 


1.0 


10 


17.257 


17.508 


17.4228 


-0.009 


0.0049 


1.0 


40 


104.317 


105.888 


105.360 


-0.009 


0.0050 
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Once again, minimizing -En (7) with respect to 7 leads to the equation, 

73-^^-2A(^^j^0. (20) 

Clearly, for small A and small n, the solution to this equation is 7 « m. In other words, 7 is what it would be 
for the usual perturbative expansion. Note, however, that no matter how small A is, the term Xn(n — 1) eventually 
dominates the energy denominator and shuts off the A'^l-growth of the perturbation expansion. Also, for large enough 
n, it is clear that Eq. ^]is eventually dominated by the term proportional to A, and for those states the energy is 
proportional to X^^^n'^^^. 

As in the previous case, numerical studies show that independent of A, m and n the second order perturbation 
theory determined by the variational calculation gives the eigenvalues of H to an accuracy of better than one percent. 



IV. EXAMPLE 3: THE DOUBLE WELL 



It should be clear from the previous discussion that adaptive perturbation theory can be applied to any Hamiltonian 
of the form 

^ = y + ^-^ + Ax^ (21) 

The problem becomes more interesting when the mass term is chosen to be negative and the potential energy develops 
two distinct minima. (This is interesting in the case of a field theory because it is associated with the sponataneous 
breaking of the discrete symmetry — )■ — 0.) I will now show how to apply adaptive perturbation theory to this 
problem, since obtaining accurate answers requires generalizing the procedure in a way that will be very important 
in the next paper in this series. 

Up to an irrelevant constant, the most general negative mass version of the anharmonic oscillator can be written 
as: 

H=lp' + lx{x'-f)\ (22) 

Clearly, for a non- vanishing value of f^, this potential has two minima located a,t x — ±/. Thus, we would expect 
that the best gaussian approximation to the ground state of this system can't be a gaussian centered at the origin; 
classical intuition would imply that it is a gaussian centered about another point, x = c. In other words, if |0^) is a 
state centered at a; = 0, it is better to adopt a trial state of the form 

\c,j)=e-'^P\0,). (23) 

Since 

e'^Pxe-'^P = x + c, (24) 

it follows that computing the expectation value of the Hamiltonian, Eq. 1221 in the state specified in Eq. [531 is the 
same as computing the expectation value of the Hamiltonian obtained by replacing the operator a; by x + c, in the 
state |0-y). In the language of the previous sections, this is equivalent to introducing the annihilation and creation 
operators A\ and as follows: 




With this definition, the normal ordered form of the Hamiltonian, Eq. I22L becomes 
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The expectation value of this Hamiltonian in the state |0^) is 

which should be minimized with respect to both 7 and c in order to define the starting point of the adaptive 
perturbation theory computation. 

Fig- n shows £(7,c) for a particular choice of A and and it exhibits what is a generic feature of any such plot: 
namely that there are three minima. To get a better feeling for what the three dimensional plot is showing us, it is 
convenient to hold c fixed and solve for the value of 7 that minimizes £"(7, c), which I will denote by 7(c), and then 
plot £"(7(0), c) for various values of A and Three such plots are shown in Fig. |21 The first plot is for a value of 
/ which is large enough so that the lowest energy is obtained for a gaussian shifted either to the right or left by the 
amount c = ±Cmin- While there is a local minimum at c = it has a higher energy than the shifted states. Things 
change as one lowers the value of /. Thus, the second plot shows that for lower / the shifted wavefunctions and the 
one centered at zero are essentially degenerate in energy. For a slightly smaller value of /, things reverse and the 
unshifted wavefunction has a lower energy than the shifted ones. 



E(^i,c) 




FIG. 1; Effective potential as a function of 7 and c showing three local minima. 

While there would seem to be nothing wrong with the situation shown in Fig.|21 it is problematic when one applies 
the same sort of analysis to negative mass 0** field theory in 1 + 1-dimensions. In this case, if one were to add a term 
like J0 to the Hamiltonian, this result would imply the existence of a first order phase transition when J reached 
some finite value. At this point the expectation value of in the ground-state would jump discontinuously from zero 
to a non-zero value. It has been rigorously shown that such a first order phase transition at a non-vanishing value of 
J cannot occur 

For this reason it is important to understand the origin of the secondary minimum at c = and how to avoid it 
ever becoming lower in energy than the flanking minima at c = ±c,nin- Fortunately the origin of the problem is easy 
to understand and fixing it is not difficult. The key is that the Hamiltonian in Eg. 1261 contains a term linear in c and 
linear in Aj^ -I- A,^ . Taking the effect of this term into account removes - or at least minimizes - the importance of, 
the secondary minimum at c = 0. In order to allow this term to play a role in the variational calculation, we have to 
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FIG. 2: Effective potential as a function of c, showing two degenerate local minima at c = icmin and one minimum at c = 0. 
The three plots are for A = 1 and different values of / and show how the minimum at c = changes from lying higher than the 
minima at c = ±Cmin to being the global minimum. If this were field theory this would be characteristic of a first order phase 
transition. 



adopt a slightly more general form of the variational wavefunction; namely, we try a state of the form 

■0(7, c, a) = e*''^ (cos(a) \0j) + sin(a) |1^)) . 



(28) 



The easiest way to minimize the expectation value of the Hamiltonian, Eq. 1261 in this state is to fix the values of 7 
and c and solve for the value of a which minimizes this expression. Then substitute this value, a(7, c), and minimize 
the resulting equation respect to 7 and c. However, minimizing the expectation value of the Hamiltonian when 7 and 
c are held fixed is equivalent to diagonalizing the 2 x 2-matrix 



^^2(0,7) 



/ (0^|il(7,c)|0^) (0^|il(7,c)|l^) 



V(l^|i/(7,c)|0^) (l^|i/(7,c)|l^) 
so it is simple to show that the lowest eigenvalue of this matrix is £(7, c); i.e., 



(0^|i/(7,c)|0^) + (l^|i/(7,c)|l^> 



;i^|i?(7,c)|l^)-(0^|iJ(7,c)|0^) 



(0^|i7(7,c)|l^)^ 



The eigenstate corresponding to this eigenvalue has the form 

1 



\Mc)) 



0. - 



Ac) 



(29) 



(30) 



(31) 
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where 

_ v/((0^|g(7, c)\0,) - {l,\H{j, c)\^))^ + 4 ((0^|g(7, + ((0^|g(7, c)|0^) - (l^|g(7, c)\l,)) 

2(0,|i/(7,c)|l,) ■ ^'2^ 

Note that x(c) is an odd function of c. The next step is to minimize £{j,c) with respect to 7 and c. Fig. |31 shows 
what a typical plot looks like for finite values of . For some values of A and / one can discern a secondary minimum 
near c = 0, but it always lies higher than the minima corresponding to non-zero values of c. 

It may seem surprising that although Cmin ^ as / ^ 0, it doesn't vanish at / = 0. Neither does the expectation 
value of the operator x, although it is much, much smaller than Cmin- This is because the energy is improved by 
moving the gaussian slightly off x — and widening the symmetric wavefunction, admixing a small amount of the 
N = 2 state. Note that if we rewrite the Hamiltonian, introducing a factor V as follows, 

H=±p' + ^X{x'~f)\ (33) 

and let V ^ 00, then Cmin ^ as / ^ 0. This happens because V acts as a classical mass for the oscillator and one 
expects that for infinite mass, the wavefunction of the oscillator should become a narrower and narrower and lie at 
the classical minimum of the potential. 

Returning to the case V — 1 and non-vanishing /, the fact that the two variational wavefunctions I7, -l-Cmin) and 
1 7, —Cmin) have the same energy means that no matter how small the overlap (7, — c,iii,i|i/ 17, -|-Cmin) they will mix 
maximally. Thus it is possible to do better by considering the symmetric and anti-symmetric wavefunctions 

= e'^P (cos(a) 10-,) + sin(a) |1-,)) + e-'^P (cos(a) |0-,) - sin(a) |1-,)) 
IV^odd) = e'^P (cos(a) |0^) + sin(a) |1^)) - e^'^P (cos(a) |0^) - sin(a) |1^)) , (34) 

and computing the expectation value of the Hamiltonian in each of these states. The formulas needed to carry out 
this computation appear in Appendix B. 
The most important results are: 

1. For large the energy splitting between the even and odd states is, as expected, exponentially small; 

2. This splitting grows as tends towards zero and, in the limit of = 0, becomes of order unity. 

These results say that while for large / it is appropriate to think of the splitting as due to tunneling between the 
degenerate minima, this is not the case when f 0. This is true even though c,nin stays finite. The reason c,nin stays 
finite in the limit / is that a non-vanishing value of Cmin produces a state that is an admixture of n = 2 and 
higher n states with the n = state. In fact, due to this admixing this calculation produces energies for both the 
ground state and first excited state of the pure anharmonic oscillator that are equal in accuracy to our earlier second 
order adaptive perturbation theory calculation. 
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A. Larger n 

Things change for larger values of n^, even for / ^ 0. Generically, what happens is that as grows, Cmin tends 
to zero, but doesn't quite get there. Thus as for the case = 0, there are still two degenerate minima corresponding 
to equal and opposite values of Cmin f^nd so one can lower the energy by forming the states 

IV'ovcn) = e'-'^P (cos(q;) \n^) + sin(a) \n + 1^)) + e'^^^ (cos(a) \n^) - sm{a) \n + 1^)) (35) 

and 

IV'odd) = e'^P (cos(a) \n^) + sin(a) \n + 1-,)) - e'^^P (cos(a) |iV^) - sin(a) \n + 1-,)) . (36) 

The result of such a computation is to show that as grows the splitting between these states grows. This follows 
directly from a computation of the overlap 

even ), (37) 

which, using the formulas in Appendix B, can be shown to grow as a function of n^. Eventually, no matter what 
value we assign to / (so long as it is finite) there is a value of n-y for which the splitting between the states |'0oven) and 
iV'odd) becomes of order unity. This is the point at which it makes no sense to talk about tunneling between states 
confined to one or the other side of the potential barrier. 

V. TUNNELING BETWEEN ASYMMETRIC MINIMA 

In the preceding section I showed how to use adaptive perturbation theory to compute the tunneling between 
degenerate minima in a symmetric potential. Clearly, while the adaptive perturbation theory discussion is a rather 
simple way to get accurate results, it is by no means the only method that can be used. However, if one wishes to 
consider the tunneling in a very asymmetric potential (see Fig. ^ then the familiar techniques are much more difficult 
to apply. In particular, consider the problem of computing the time evolution of an initial state corresponding to 
the lowest state concentrated in the left hand well. I know of no simple techniques which allow one to compute this 
time evolution accurately. In this section I show how to use the techniques of adaptive perturbation theory, plus one 
additional trick, to obtain a simple picture of how things evolve in time. 




FIG. 4: Asymmetric potential well. 

Prom the preceding discussion, it should be clear that given a specific form for the potential shown in Fig. ^ we 
can use the adaptive trick of rewriting the problem in terms of shifted Fock states. This allows us to approximately 
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compute the lowest lying state in the left well and a finite but sufficiently large number of n-particle states in the 
right well. These computations, augmented by perturbation theory, take care of computing the effects which have 
nothing to do with the tunneling. 

Tunneling contributions are obtained by computing the transition matrix elements between the lowest lying state 
shifted to the left of the origin and all the states corresponding to a shift to the right of the origin. These matrix 
elements are computed as in the case of the symmetric potential discussed in the preceding section. Assuming that 
this step has been carried out to the desired degree of accuracy, and that the minima are widely separated so that 
the transition matrix elements are small, the next step is to compute the time evolution of the state on the left when 
the tunneling is taken into account. 



A. The Simplified Problem 



At this point we have effectively mapped the original problem into a simpler one; namely, that of a Hamiltonian 
H which can be decomposed into H — Hq + V, where Hq has a discrete set of eigenstates {ipj) with energies Ej, 
1 < j < N, (representing the levels in the right hand well) and a single state IV'o), with energy Eq (representing 
the lowest lying state of the left hand well). Moreover, given the potential in Fig.^ it will be assumed that Eq lies 
somewhere in the range of the Ej, i.e., Ei < Eq < En, and that it is not degenerate with any of the Ej^s. This is an 
assumption I will relax at a later point in this discussion. Note that in this form the problem is essentially equivalent 
to the Wigner-WeisskopffQ problem of a discrete state decaying to a continuum, except that now it is a discrete state 
decaying to a large number of discrete states. Since I have assumed that adaptive perturbation theory was used to 
determine the Ej's and Eq, the perturbation V will be obtained by computing transition matrix elements between the 
state {ipo) and the states \ipj). In other words the only non-vanishing matrix elements of V will be Voj — (i/'ol^lV'j) 
for l<j<N. 

I will now discuss how to compute the probability of finding the system in the state lipo) or lipj) at any time 
t, assuming it is in the state \tpo) at time t — 0. This is, of course, equivalent to computing (V'o|e~**^|^o) or 



B. The Resolvent Operator 

For reasons which will become clear, the best way to calculate (■0o|e**'^|'0o) for i > 0, is in terms of the resolvent 
operator 1/{H — z). This is because the matrix elements of the resolvent operator are directly related to the Laplace 
transform of the function in question. Explicitly, consider the Laplace transform 

V'oo(^) = r dte-^' {tPj\e-''"\iJo), (38) 
Jo 

where z = c + ico for fixed c > and — cxi < a; < oo- Because, in our simplified problem, the eigenvalues of H are real 
and finite in number, this function is obviously well defined for z in the right-half complex z-plane, and for t > the 
inverse transform is given by the contour integral 

-1 pc+ca 

{^j\e'^*"\^b,)^ i:,o{z)e''dz. (39) 

^TTl Jc-ioo 

The function is defined in the full complex plane by analytic continuation. 
Carrying out the integration in Eg. I38lvields 

V'oo(^) - (40) 
H — iz 

Since the resolvent operator has poles onlyl at the discrete eigenvalues of the operator H, it is convenient to substitute 
z — > iz and redefine the contom in Eg. 1391 to obtain 

(V|e'*^|V) = - / dz(V|-— |V')e-^*^ (41) 

where we have assumed t > 0. Obviously, the same thing can be done for two different states; i.e., 

— oo+ic -1 

(x|e'*^|V^) = - / dz (xl-— |^)e-^*^ (42) 
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The reason for focusing on the resolvent operator is that it will be easy to compute the relevant matrix elements of 
the resolvent operator exactly. This is because the resolvent operator satisfies the equation 

-^^^ '—V^ (43) 

and therefore the matrix element of interest satisfies the corresponding equation 

(V-Itt^I^) = (V^Itf^IV') - Wj7^V-^\i,), (44) 

H — Z IJQ — Z Hq — Z 11 — z 

which can be solved by iteration to obtain the formal solution 

1 1 °° 1 

(V-Iti— IV'> = (V'|-5-— E(-^i:r^)"IV')- (45) 

tl — z tlQ — z ^ rlo — z 

n— 

For our problem, this sum can be evaluated in closed form because the state l-i/^o) is an eigenstate of Hq with 
eigenvalue Eq and the only non-vanishing matrix elements of V are of the form V^j ~ (V'ol^l'^n)- From this it follows 
that the sum in Ea. 1451 reduces to a geometric series, which is easily summed to give 

This can be cast into a more useful form by multiplying both the numerator and denominator of the expression 
by Jlj(^j ~ ^) ^'^'i rewriting it as a rational polynomial. Then the poles of the resolvent operator correspond to the 
zeros of the denominator and so we can rewrite the resulting expression as 

1 1 nf-i(^^7-^) 

Note that here the eigenvalues of H have been denoted as Ej. It will be convenient in what follows to rewrite the 
Ej^s in terms of their difference from the corresponding eigenvalues of Hq\ i.e., 

E,^E,+5,. (48) 

It is the quantities Sj which need to be determined. 

Eq. E3is exact, independent of the size of V. However, since the matrix elements of V are assumed to be quite 
small, it is useful to rewrite it in a way that allows the Sj^s to be computed perturbatively. To do this, first rewrite 
the second term in Eq. 03 as 

N N N 

l[{E,-z) = l[{E,-z)-\Voifl[(E,-z) 

N p-1 N 

- \Vo2\'{E,~z)l[{E,-z)- ...-iVo.l'lliE.-z) l[{Ei-z) 

3=2 j=l l=p+l 

- ...-\Von\'1[{E3-z). (49) 

i=i 

From this we see that |VbjP ^ 1 implies Sj w |VbjP- Thus, for example, setting j = and choosing z = Eq + Sq, the 
equation becomes 

N N N 

= -So HiE, - Eo) - \Voif n -1^021' {El ~ Eo) J] 

3 = 1 3=2 j=3 

p-1 N N-1 

- 1 n(^^- - n - - . . . - \von\^ n (^^- - ^o), (5o) 
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or in other words, 

j=l J 

The other energy shifts are obtained in a similar manner. For example, if we set z = E\ — E\ ^ 6\ then, to lowest 
order in b\ or | Vg j | ^ , the equation becomes 

JV JV 

This is true because the other terms take the form 

JV m-l JV 

- i^02p (-<5i) n(^j- - - . . . - iFo™p n (^.- - ^i) n (^j- - ^i)' (53) 

which are all of second order in the perturbation. Thus, finally, we have 

^1 = (54) 



Similarly, we for all other values of j we obtain 



- (55) 



(Note, I have implicitly assumed that the spacing between the levels Ej is larger than the typical value of |VojP 
and that Eq is not degenerate with any of the -Ej's. The second restriction is easy to remove, since if Eq k, Ej for 
some j, then computing Sq and bj just requires solving a quadratic equation for its two roots. The other (5j's remain 
unchanged. If however the spectrum of the EjS becomes dense it is necessary to treat the problem of computing the 
levels which lie within a distance of order |VbjP of Eq more carefully.) 

These results are half of what is needed to reconstruct the full time dependence of ('!/'o|e'*^|'0o)- It follows from 
Eq. 1391 that since the matrix element of the resolvent operator is a sum of simple poles, we can get the full time 
dependence of the matrix element by closing the contour. This would give the time dependent behavior as a sum of 
exponentials multiplied by the residue of the corresponding pole. Fortunately, given the previous discussion, it is easy 
to compute the residues of the poles. To compute these residues we begin by rewriting Eg. 1471 as 

^ " ^ Eq-z-Y^j ^e73I 11^=0 \E^ - z) n.=o yEs-z) 

and then expand the result to first order in 5j . This yields the equation 

Now observe that since Eq. |23 explicitly shows that there is a pole at every value z = Ej it becomes a simple matter 
to compute the residues to the same order in the perturbation. Letting z Eq yields the result 



and, similarly, letting z — > Ej implies 
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Combining these results, we see that the amphtude for finding the particle in the state |'0o)j at any time t, is given 

by 

N 

(V'o|e^^*|^o)=E-^J-e"'''"^*- (60) 
It follows from evaluating Eq. (HOI at t = that 

QED. 

To get the amplitude for making a transition to state iV'm): we have to compute 

(V'^Itt^I^o). (62) 

However, since V only causes a transition fr'om ji/'m) to IV'o); it follows fr'om the integral equation 

1 _ 1 1 1 

H - z ^ Ho- z ^ Ho- z ^ H - z ^^^^ 

that 

i^m\-n—\^0) = -J^ 'at • (64) 

H-Z hrn~Z \lll=oiEl - Z) J 

Canceling the factors of {Ej — z) wc obtain 

(V'mlTJ \^o) = V,nO— ^ • (65) 



Expanding this expression to first order in the 5's we obtain 



Em — Eo ^ {Em — Eo) {Ej — Eq) ;J^^ {Ern — Eq) {El — Eo) 

TA m— 1 r. ^ X 

_ VmO \ ^ \ - 01 



E, 



(5 

fj = 7^ ^TTF where Oand TO. (66) 

It is important to check that this solution gives the correct boundary condition that at i = 0; i.e., (?/'m| e'^* |V'o) — 0- 
Clearly, this is the same as requiring that the sum of the residues vanish. Applying the preceding result we see 
that 

5, / 1 \ \ , ^ ( \ 1 



. — -E^O ~ -^0 -^7 ^ / - — Eq \Ej — Eq Ej ~ Ej^ 

m— 1 r- N 

= y ^ ^ + y ^ ^ ^, (67) 

^ {E, - Eo) {E, - Em) ^^^^ {El - Eo) {Ei -Em)' 

and therefore, 

m-l N 

ro + f„ + E fj + ^' = 0, (68) 
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which is the desired result. 

Given these formulas, we can write the time-dependent state |V'o(i)) as 

N 

\Mt)) = e^'"\i^o) = E (V'm|e^*''lV'o)|V'm). (69) 

m— 

It is easy to see that in the case where the eigenvalues of Hq are widely spaced on the scale of |Vbjp, the largest 
residues, and therefore the largest contribution to \ipQ{t)), comes from \ipo) and the state |'0m) for which Em — Eq 
takes the smallest value. However taking the limit in which the levels in the right well become dense, as would be the 
case if the right hand minimum is very flat and very wide, cause things to change. While the formula Eg. 1461 remains 
exact the sum over levels is best approximated by an integral over a continuous variable E; i.e., 

Since in this case Eq lies in the region of integration, we see immediately from the fact that the integral develops an 
imaginary part for Emm ^ z < Emax, that the resolvent operator no longer has a pole at real z ^ Eq. Instead, the 
zero of the denominator, and thus the pole of the resolvent operator, moves off to a complex value of z on the second 
sheet of the function. In order to compute the time dependence of the state one has to deform the contour for the 
inverse transform and pick up this pole. This will give the usual exponential decay of the state |V'o) plus a background 
integral that does not decay exponentially in t for large t. 



VI. A LOOK AHEAD TO FIELD THEORY 



Before concluding, it is worth taking a quick look ahead to the case of an interacting scalar field theory, to see how the 
techniques of adaptive perturbation theory can be used to non-perturbatively derive the structure of renormalization. 
In this section I focus on a simple calculation that gives a non-perturbative approximation to mass renormalization 
and I will leave the discussion of wavefunction and coupling constant renormalization to the next paper. 

To simplify the discussion I will limit myself to the case of a self-interacting scalar field, with a Hamiltonian of the 
form 



(71) 



The spatial dimension n is arbitrary. In what follows I assume that everything is defined in a finite volume, V , so 
that momenta are discrete. This isn't necessary, but it will make all of the manipulations which follow well-defined. 
Given this assumption, the Fourier transforms of the field and its conjugate momentum are defined to be 



m = / d"xe-^'=-^>(x) 



n^(fc) = / d"xe-^^-%(a;). 
vV . 



(72) 



In terms of these operators the Hamiltonian can be rewritten as 



H^(-fc)H^(fc) {k-k + m^) 



m<t>{-k) 



X 1 

6V 



E -^(^l) <^(fe)^(fc3) 0(fc4)<5(E ^0- 



i=l 



(73) 



At this point I set up the adaptive perturbation theory calculation by defining 7(fc)-dependent annihilation and 
creations operators as follows: 



m 



{A{-k)^ - A{k)) ; H(fc) 



(Ai-k)^ + A{k)) 



2 7(fc) 



(74) 



and define the associated vacuum state for this choice of 7(fc)'s, 



\vac) 



k 



(75) 
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by the condition that it be annihilated by all the A{kys. 

Adaptive perturbation theory around this vacuum state is now determined by minimizing the vacuum expectation 
value of the Hamiltonian with respect to all of the 7(A;)'s. It follows directly from these definitions that the function 
to be minimized is 



{vac\ H \ vac) — 



7(fc) {k ■ k + m^) 
4 47(fc) 



A 



(76) 



If the range of the momenta appearing in these sums are unrestricted, these expressions diverge, as happens in all field 
theories. The usual way this problem is dealt with in the context of ordinary perturbation theory is to regulate the 
integrals and then add counterterms to the Lagrangian to cancel the divergences. Since I wish to discuss this theory 
non-perturbatively, I will adopt a different strategy. I will render the theory well defined by restricting the operators 
(j){k) and Tl^[k) to be finite in nmubcr. This can be accomplished in a variety of ways. One is to put the theory on 
a lattice by defining the momentum components to be fc^(j) = 27r j'/L^a, where is the size of the lattice in the /i 
direction, a is the lattice spacing and j is an integer running from —L^/ a < j < L^/a. Another way to accomplish 

the same thing is to keep only those operators 4>{k) T\ci,{k) for which k ■ k < A. Since nothing I will say from this 
point on requires that I specify a particular scheme, I avoid making any detailed assumptions. 
Minimizing Ea. 1761 with respect to each j{k) yields 



1 
4 



{k ■ k + rn?) 
47(fc)2 



A 



27(fc)2 



1 



which can be rearranged to give 



k-k 



2A 



-y ^ 



V ^ 7(fc') 



(77) 



(78) 



Setting fc = we obtain 



7(0)^ = m2 + 2A 
which can be substituted into Eq. [7H|to give 



-E^ 

V ^ 7(fc') 



7(fc)2 = fc.fc + 7(0)2. 

Substituting this expression into the equation for 7(0) yields the non-perturbative equation 



7(0)' = m2 + 2Al5: 



(79) 



(80) 



(81) 



fc, -v/fc'.fc' + 7(0)2 



which is reminiscent of the Nambu Jona-Lasinio equation. This equation is guaranteed to have a solution since the 
formula for the vacuum energy clearly diverges both at 7(0) = and 7(0) — 00 and is finite in between. At this point 
we can take the infinite volume limit and replace the sum over fc's by an integral up to some cut-off A. In this case 
the equation to be solved becomes 



7(0)= 



2A- 



2P 
(2^ 



kP-^dk 



1 



k • fc + 7(0)2 



(82) 



The simplest way to understand the general properties of the solution to this equation is to graph the left and right 
hand sides as a function of 7(0) and see where they cross. Clearly, the graph of the left hand side of Eq. |221is a 
parabola that intersects the vertical axis at —rri^. The graph of the right hand side of the equation starts at some 
large (or in l+l-dimension infinite) value and then drops monotonically to zero as 7(0) — > 00. Generically, since for 
finite 7(0) <C A, the integral is proportional to A^"^ for p > 1, and ln(A) for p = 1, the point where the parabola 
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crosses the graph of the integral will be for large 7(0). The only way to have the graphs cross for a value of 7(0) 
that is small with respect to A is to make large and negative. In fact, one should really think of Eq. |S2not as an 
equation for 7(0), but rather as an equation for m^; i.e., rewrite it as 



and then arbitrarily pick a value of 7(0) and use this equation to determine for what value of the chosen value of 
7(0) will minimize the ground state energy density. This is nothing but a non-perturbative way of determining the 
leading mass renormalization counter term. 

While this discussion is amusing, it does nothing to show what one has to do to similarly understand wavefunction 
and coupling constant renormalization in a non-perturbative context. That discussion will be the subject of the second 
paper in this series. 



It was my purpose in this paper to introduce the notion of adaptive perturbation theory as a way of perturbatively 
computing quantities that heretofore couldn't be computed in any simple way. I focused on problems in ordinary 
quantum mechanics such as the pure anharmonic oscillator and the problem of computing the eigenstates of generic 
double-well problems in order to exhibit the versatility of the technique. Finally, I discussed the simplest way to begin 
to apply these techniques to a field theory in order to obtain a better behaved perturbation expansion independent 
of the size of the interaction term. In addition, I tried to indicate how the same approach would allow one to non- 
perturbatively approach the computation of renormalized quantities. In the next paper I will show how to extend 
these techniques to non-perturbatively capture the full structure of renormalization. I will also show how something 
related to the running coupling constant appears and explain why this approach essentially forces one to discuss the 
computation of scattering amplitudes in a way that is intimately related to the approach used in the parton model. 

Obviously there are many items which I have not touched upon and which require closer analysis. One particularly 
interesting unsolved problem, which is relevant to the analysis of a compact gauge theory, has to do with a problem 
in which the potential has an infinite number of degenerate minima. When these minima are widely separated by 
barriers of sufficient height then one can do a good job simply by using the adaptive technique to compute states 
localized in any one minimum and then summing over appropriate superpositions of these states to obtain Bloch wave 
function. However when there is significant tunneling between nearby minima one has to better understand, as we 
did in the simple tunneling problem, how to include single-well excited states in the trial wavefunction so as to avoid 
double counting in the perturbative expansion. 

The key point I would like the reader to take away from this discussion is that, in an interacting theory, the 
perturbation theory one uses to compute an interesting quantity must be adapted to the problem at hand. Although 
a variational computation based on energy considerations is the tool used to define the appropriate parameters to be 
used when dividing the Hamiltonian up into an unperturbed part and and an interaction term, it is the perturbation 
theory based upon this division which is crucial to computing quantities other than the energy of a state. This 
observation takes us beyond the more usual variational methods which cannot be usefully applied to compute matrix 
elements of general operators, or scattering amplitudes in a field theory. 



I would like to thank Carl Bender for useful communications and in particular for pointing out to me the existence 
of the papers by Halliday and Suranyi|2(|. 



This appendix is devoted to reminding the reader how to derive normal ordering identities that are useful for 
carrying out the computations I have already described, as well as the computations in Appendix B. 




(83) 



VII. REPRISE 



VIII. ACKNOWLEDGEMENTS 



APPENDIX A. SOME OPERATOR IDENTITES 



Theorem: 



^aA\ ^aA., 



(84) 
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Proof: Assume 

and differentiate with respect to a to obtain 

(4 + A^) e^f^U-^T) ^ [ + ^5(a)4 + -^h{a) e^^^^^^ e'*^")^- g-^^")^^! e^(") e^^")^' e'^'")^"'. (86) 

' \_da da da J 

Successive differentiation of the term e^*-"^'^^' e''^"^"^^ e"^'-"-'"^^ with respect to a gives the identity 

gSMAt e-«W^^ = Yl ^^K' [4' [■ ■ • ' K' A]]] ■ ■ •]"■ (87) 

n 

Observing that [v43y,v4^] = — 1 ahows us to rewrite this as 

gS(a)At gh(a)A^ g-ff(a)At ^ _ g(^^y (gg) 

Substitute this into Ea.lSBIand cancel the factor of e"^'^^"'"'^^-' against the term e^*^"' g9{a)A\ ^h(a)A^ ^j^g right hand 
side of the equation to obtain 

(4 + A) = ^/(") + ^SH^ + ^^(") (A - • (89) 
Equating coefficients of Al^, A^ and the unit operator yields the following differential equations: 

^g{a) = l ; ^h{a) = l ; ±f{a)=g{a)^h{a) (90) 
da da da da 

The initial conditions /(O) = 0,5(0) = and h{Q) = follow from the fact that e" the identity operator for 

a = 0. Given these boundary conditions, it follows directly that 

a^ 

g{a)^a ; h{a) ^ a ; f{a) = —, (91) 

which is what we wished to prove. 

A similar calculation proves the following result: 
Theorem: 

This identity is used to compute, among other things, the overlap (0-y|e~*^^|07) as follows: 

{Q^\e-'^P\0-,) = (0^|e^\/*(<"^^)|0^) =e-^'^/4(0^|e=V¥4e"VI^-|0^) = e-"'T/^ (93) 

since, by assumption, yl^|0-y) = 

I mentioned that there was a simple expression for the state |0^) in terms of the states built upon 10-^'). The explicit 
form of this relation is the following: 

Theorem: 

|0^) = VT^teS^e-^ '""'^('')^''|0y), (94) 

where tanh(/3) = (7 — 7')/ (7 + 7')- 

Proof: Given that we introduced that 7-dependent annihilation and creation operators as 



A\+A^ [2 

^ = — ; p 



; p = ^^^{Al-A,) (95) 
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it follows that for two different values of 7 



2 vtt' 



By definition A^|0-y) =0. Thus we have 



Now, consider a state of the form 



e 1' |0-y'). 



Requiring that annihilates this state is equivalent to the statement 

1(7 + 70 . . A(At )^,^ ^ , 1(7-Y) .t A(At ) 



2 



2 



4^e^^^>^ |0y) = 0. 



AA+ 



Moving the factor e t^' to the left of each term we obtain 



e T 



2 1/7 7' 2 ■v/77' 



|o;) = 0. 



(96) 



(97) 



(98) 



(99) 



(100) 



Once again, expanding the term e -v' A^i e -v' in multiple commutators we obtain the results that will 
annihilate the state if and only if 



or, in other words, if 



QED. 



(101) 



(102) 



Finally, the normalization factor come from the identity 



e'^''"|0y) = 5:^V^|2ny) 



A" 



n=0 



The norm squared of this state is given by the sum 



Vl - 4A2 ' 



so the normalized state is 



|0^) = (1 -tanh(/J)2)V4e**""^^''^^v'|0y). 
It is worth noting that another way to derive this result is to observe that 

is a unitary operator that satisfies the relation 

-|(At'-A2) ^ -i tanh(/3)At^ - l„(cosh(/3)) (A+A+i) i tanh(/3)At' 



(103) 



(104) 



(105) 



(106) 



(107) 



This is proven in much the same way as are the other normal ordering relations. Putting back the appropriate 7' 
dependence in this equation and applying the operator to the state |0-^') gives the desired result. 
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APPENDIX B. DETAILS OF DOUBLE WELL CALCULATION 



As we already noted, the two degenerate variational states for non-vanishing are defined as: 

|x(c,7))=e-^^P|^^(c)) (108) 

and 

|x(-c,7))=e*^P|V^(-c)) (109) 

(assuming c > throughout). 

Once one has found the best choice of variational parameters, still lower energy states are obtained by forming the 
even and odd parity combinations: 



|^oven(7)) = T7^ ( lx(c, 7)) + lx(-C, 7)) ) (HO) 
'even 



and 

|V'odd(7)) = TT^ ( lx(c, 7)) - lx(-c, 7)) ) , (111) 

JVodd 

where N^vcn and A^o^d a-re normalization factors. 

To evaluate the normalization factors Ncvcn and A'odd we need to compute the scalar product (x(c, 7)|x(— c, 7)) as 
well as overlap elements {x{c,j)\H\x{—c,"f)). Since the diagonal scalar products arc obviously unity, the diagonal 
expectation values of the Hamiltonian are simply the values previously computed by minimizing £(0,7). 

Expressing p as a sum of annihilation and creation operators and using the identity for computing the normal 
ordered exponential, Eq. 1921 allows us to rewrite the overlap as 

(x(-c,7)lx(c,7)) = {^J-c{j)\e-^''P\Ml)) = e-^'^(V-c(7)|e=^^' e-^^^^\Mj)) = e"^'^(0_c(7)l0c(7)), (112) 
where the states (0c(7)| and |(/)_c(7)) are defined in terms of the states 

|Vc(7)) = ,^ \ IMl)) - lV'i(7)) (113) 

v/1 + 2:(c)2 ^l + x{c)^ 

and 

Note that in this formula x{c) is given by the formula evaluated for positive c, and the states |0c(7)) and |0_c(7)) 
are given by 

|0c(7)) = e-^^^-' \Ml)) (115) 

and 

\ct>-ch))^e^^^^H^ch)). (116) 

Since the state lipdl)) contains only n = and n — 1 states, only the first two terms in the expansion of the 
exponential act. Thus it follows that 

I0c(7)) = e~^^^-^\Ml)) = ( ^^p^^^A 1^0(7)) - -jM=\Ml)) (117) 
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and 



so the overlap is 



(X(-C,7)lx(c,7)) =e-''{0_.(7)|*,(7)) =e-"-^ii^i^5*al^. (119) 
Thus, the normalization for the state |'^even(7)) is 

A^even = V + {x{-C, iMc, j))) . (120) 

Similarly 

iVodd = V2(l-(x(-c,7)lx(c,7)))- (121) 

Finally, to compute the formula for the expectation value of H in the state |'^i'even(c, 7)) or |V'odd(c, 7)), we need the 
following identities: 

e-'"Pfl'(p,a;)e'"f = H{p,x-c) 

e^^P H{p,x)e-^''P = H{p,x + c) 

e-'^PH{p,x)e-'^P = H(p,x-c)e-^"'P = H{p,x-c)e^^^^'re-^^^-^ 

= e-'='^e'=^'^^7?(p-ic7,a;)e-"^^'^^ 

e"'PH{p,x)e"'P = iJ(p, X + c) e^^cf = e-"''' H(p, .T + c) 6^=^^' 6^=^^^ 

= e-'-'-'e-''^^^' i^(p + ^c7,a;)e"^'^^. (122) 

We then need to evaluate these expressions between the states I'fpdl)) and \tjj-c{'y))- Clearly, the first two terms are 
equal to £{X, P,c, 7)inin, which is the minimum value of £{X, p, c, 7) obtained by minimizing with respect to c and 
7, since the minima are degenerate. 

Given these identities, it is straightforward to evaluate 

(V'even|i?|V'even) = "TT^ [{x{c, j)\H {p, x)\x{c, -f)) + {x{-cn)\H{p, x)\x{~C,-f)) 
even 

+ {x{-c,l)\Hip,x)\x{c,j)) + {x{c,^)\H{p,x)\xi-c,^))]. (123) 
Clearly, the first two terms evaluate to the function £{X, p, c, 7) and the remaining terms are 

{x{-c,j)\H{p,x)\x{cn)) = e-'='^((l + c^A^)^_,{'y)\H{p-in,x)\{l - c^A^)Ml)) (124) 

and 

(X(c, j)\H{p, x)\xi-c, 7)) = e-^'^((l - c^A^)Ml)\\Hip + icj, x)\{l + c^/2^A^)V-c(7))- (125) 
Explicit substitution of the shifted definition of p in iJ yields 

H{p — ic^, x) = H{p, x) — icjp — —■-f (126) 



and 



.,2 



H{p + icj,x) = H{p,x) + ic^p — ^7- (127) 
To relate the two cross-terms, it is helpful to observe that 

(x(-c,7)|0|x(c,7)) = (X(c,7)|0^lx(-c,7)>*- (128) 
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Since ip is anti-hermitian, we get a minus sign for that term and plus signs for all the others. Note that the expectation 
value of H{p, x) that appears in the cross terms is obtained from Ea. l29l sctting c — and evaluating the result between 
the states (0_c(7)| and \(f>c{'y))- 

The result for the expectation value of the Hamiltonian in the variational states IV'even) and jV'odd) is: 



(V'ovenl^^lV'cvon) = 



l + a;(c)2 



il + y^cx{c)) (t + — + ^-^ 



+ 2x{c)il + ^2jcx{c))cj\jj-x{cy 



and 



5A 37 XP , A/4 cV 
872 4 



27 



6 



(129) 



(V'oddl^^lV'odd) 



iV2 



^(A,/^c,7) 



-c 7 



l + a;(c)2 



+ 2x{c){l + ^cx{c))c^^l-x{cf 1^— + ^-^ 



7 A/^ A 
4 6 872 67 

5A 37 A/2 A/4 C272 



2 2 

c 7 



6 



(130) 



So long as the factor e~'^ is small, the splitting between the even and odd states is quite small. Thus the tunneling 
between the left and right states is slow. When we go to smaller values of / or larger values of n, the tunneling rate 
increases and even at / = minimizing (V'evenl^^lV'even) and (V'oddl-ff IV'odd) with respect to 7 and c leads to a non- 
vanishing value for c. The energies obtained in this way agree in accuracy with the second order adaptive perturbation 
theory calculation for the pure anharmonic oscillator, which was done without using the shifted wavefunction. As 
expected, this shows that the effect of having a non- vanishing value of c at / = is to include the effects of the second 
order computation. 
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this computation predicted a first order phase transition at a non-vanishing value of external field and Coleman pointed 
out to us that, in the case of 1 -|- 1 dimensions, this violated a theorem by Simon and Griffiths 5]. At the time we didn't 
really understand the origin of this problem. Subsequently the same calculation was published by S. J. Chang Phys. Rev. D 
12 (1071)1975, where he too pointed out that the result violated an exact theory in the case of 1 -I- 1 dimensions. Chang's 
paper was followed by a paper by Paolo Cea and Luigi Tedesco Phys. Rev. D 55 (4967), which presented a field theory 
improvement of S. J. Chang's computation based upon a diagrammatic perturbation theory which eliminated by false 
tri-critical behavior, but not at the level of the variational computation. Subsequently a systematic perturbation theory 
about the variational calculation was discussed by Jeff Rabin Nuclear Physics B 224 (308-328) 1983. So far as I know, no 
variational way of removing the tri-critical behavior was published before the analysis in this paper. 



